Higher order QED in high mass e + e pairs production at RHIC 
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Lowest order and higher order QED calculations have been carried out for the RHIC high mass 
e + e~ pairs observed by PHENIX with single ZDC triggers. The lowest order QED results for 
the experimental acceptance are about two standard deviations larger than the PHENIX data. 
Corresponding higher order QED calculations are within one standard deviation of the data. 
PACS: 25.75.-q, 34.90. +q 
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A recent publication of the PHENIX collaboration at 
RHIC has presented data on ultraperipheral Au + Au 
photoproduction of the J/ip and of continuum e + e~ pairs 
in the invariant mass range of 2.0 - 2.8 GeV/c 2 [H. The 
present paper concerns itself not with the J/vb but with 
the continuum pairs. A previous measurement of e + e~ 
pairs at RHIC was carried out by the STAR colaboration 
in the somewhat lower invariant mass range 140 - 265 
MeV/c 2 0, and it was subsequently argued that these 
data exhibited evidence for higher order QED effects Q. 
This paper presents results of a higher order QED calcu- 
lation for the PHENIX data using the same methodology 
as the calculations [3] previously published for the STAR 
data. 

The PHENIX data was presented as a differential cross 
section with respect to the pair equivalent mass m e + e - 
and the pair rapidity y pa ir with the constraint of at least 
one neutron detected in one of the zero degree calorime- 
ters (ZDC): d 2 o-/dm e + e -dy pa ir(e + e~ + Xn,y pa i r = 0). 
The differential cross section at y pa ir = was taken from 
the sample \y pa ir \ < 0.35 but with no constraint on pseu- 
dorapidity ry of individual electrons and positrons. How- 
ever, since the events measured by PHENIX detector had 
both the electron and positron pseudorapidities within 
\r)\ < 0.35, the stated results without the individual 
\r)\ < 0.35 were necessarily constructed from a model cal- 
culation. The PHENIX publication states, "The fraction 
of events with \(y pa ir)\ < 0.35 and 2.0 < m e + e - < 2.8 
GeV/c 2 , where both electron and positron are within 
\r/\ < 0.35 is 1.10%. The corresponding numbers for 
2.0 < m e+e - < 2.3 GeV/c 2 and 2.3 < m e+e - < 2.8 
GeV/c 2 are 1.11% and 1.08%, respectively." Since we 
wanted to carry out calculations closest to what PHENIX 
actually measured, here we calculate the cross sections 
with the electron and positron constraints |ry| < 0.35 and 
compared with the cross sections in the PHENIX publi- 
cation multiplied by 0.011, etc. (see also Ref. [4]). 

Following the method of the higher order calculations 
previously used for the STAR data, the cross sections 
here are computed from the product of the pair produc- 
tion probability P ee (b), the probability at least one of the 
ions Coulomb dissociating P lx (b), and a survival factor 
exp[—P nn (b)] to exclude events where hadronic interac- 



tions occur 



° = 2k J P lx {b)P ee (b) exp[-P nn (b)]bdb. 



(1) 



Unlike the STAR Coulomb dissociation factor P xx (b), 
which corresponds to neutrons detected in both ZDCs, 



P xx (b) = [1 - exp(-P c {b)} 



(2) 



here the corresponding Coulomb dissociation factor 
P\x{b) is the unitarized probability that requires only 
that at least one of the colliding nuclei suffers Coulomb 
dissociation 



Pix(b) = [1 - exp(-2P c (6)] 



(3) 



Pc{b) the non- unitarized probability of a single Coulomb 
excitation calculated in a phenomenological model for 
the ZDC neutrons derived from photodissociation data[5, 
6]. The non- unitarized hadronic interaction probability 
Pnn{b) is calculated in the usual Glauber manner 

Pnn{b)=<J nn J dxdy T A (x,y) T A (x - b,y)}. (4) 

where a nn is the total hadronic interaction cross section, 
52 mb at RHIC, b is the impact paramete, and the nuclear 
thickness function T A (x' , y) is the longitudinal integral of 
the nuclear density, p{r) 



T A (x',y) = / dzp(x\y,z)dz. 



(5) 



The nuclear density profile p(r — y x' 2 + y 2 + z 2 ) of 
heavy nuclei is well described with a Woods-Saxon dis- 
tribution. We use parameters determined from electron 
scattering data (R=6.38 fm a=.535 for Au). 

The perturbative and higher order QED pair produc- 
tion probabilities P ee (b) were calculated using the meth- 
ods previously described 0, 0|- I n the Pee(b) calculations 
an analytical elastic form factor was employed 



/(?) = 



{qrf 



[sin(gr) — qr cos(qr)] 



1 



1 + a 2 q 2 



(6) 



with a hard sphere radius r = 6.5 fm and Yukawa poten- 
tial range a = 0.7 fm. This form very closely reproduces 
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TABLE I: RHIC: Au + Au, 7 = 107, single ZDC triggered 
e + e~ pairs cross section (da/dm e + e - dy pair [fib/GeV]) with 
PHENIX cuts (see text). 





Data 


Starlig 


tit Perturbative Higher Order 


(GeV/c 2 ) 


[1,J 


[9] 


QED [3, 7] 


QED [3, 7] 


2.0 - 2.8 


0.95 ± 0.31 


0.99 


1.69 


1.24 


2.0 - 2.3 


1.43 ± 0.61 


1.53 


2.60 


1.90 


2.3 - 2.8 


0.65 ± 0.30 


0.66 


1.14 


0.83 



the Fourier transformation of the Au density with the 
Woods-Saxon parameters mentioned above[8|. 

Calculations in Ref.[l| performed with the computer 
program Starlight differed from the present calcula- 
tion only in the use of an equivalent photon probabil- 
ity P£. w (b) in Eq. (1) instead of the QED probabil- 
ity P ee (b). The Coulomb dissociation probability P\ x {b) 
and hadronic survival probabilities exp[— P nn (b)] were 
treated as described above identically in both calcula- 
tions. The equivalent photon probability P™ e w (b) uses 
the Weizsacker- Williams photon spectrum calculated for 
the heavy ions and the Breit- Wheeler cross section for 
(7(7 + 7 — > < ' ( ) [9j . In Starlight it is also required that 
photons originate only outside the nuclear radius R of 
each of the ions, and in the PHENIX calculations, a hard 
sphere radius R = 6.98 fm was used for this sharp cutoff. 

Results are shown in Table I. The stated statistical 
and systematic errors in the data have been combined in 
quadrature. As indicated above, calculations are for the 
measured Au + Au reaction at 7 = 107 for each of the 
Au nuclei. e + e~ pairs were accepted for \y pa ir \ < 0.35 
with a ZDC trigger indicating simultaneous dissociation 
of at least one of the gold nuclei. Individual electrons 
and positrons were also accepted only for mid-rapidity 
\r)\ < 0.35. The calculated cross section is defined for a 
pair to literally fall within the cuts: y pa ir, individual e + 
and e~ \rj\, and m e + e -, but normalized to unit y pa ir and 
m e+e -. 

As was shown in the original PHENIX publication, 
the Starlight calculations are in very good agreement 
with the data. The present QED perturbative calcula- 
tions are about two standard deviations higher than the 
data, while the higher order calculations are higher than 
the data but within one standard deviation for all three 
m e + e - ranges. 

One might reasonably ask why the peturbative QED 
calculations are so much higher than the perturbative 
Weizsacker- Williams calculations. In the first place one 
has to recognize that Weizsacker- Williams calculations 
are an appproximation to the more proper QED calcu- 
lations, in particlar neglecting the viruality of the in- 



termediate photons. However there is another physics 
difference between the two calculations. In the QED cal- 
culation there is no sharp cutoff analagous to the cutoff at 
the hard sphere nuclear radii R in the STAR equivalent 
photon calculation, but the form factor, Eq. (6) gives a 
smooth cutoff in both tranverse and longitudital photon 
momenta arising from the nuclear Woods-Saxon charge 
distribution. This smooth cutoff allows more high mo- 
mentum photons contributing to the cross section than a 
sharp cutoff does and should increase the cross section. 

A simple exploratory luninosity calculation can be car- 
ried out for the case of a uniform hard sphere nuclear 
charge distribution by allowing contributions from within 
the ions but with an effective Z multiplied by (bi/R) 3 , 
where bi is the distance of the photon source from the 
center of the nucleus. The result is that with this finite 
density cutoff the two photon luminosity in the STAR 
measured equivalent pair mass range would increase by 
44% over the sharp cutoff, and one can infer a propor- 
tional increase in the predicted perturbative Weizsacker- 
Williams pair cross sections. 

The physical question is whether equivalent photons 
arising from within one the ions when the ions are not 
overlaping should contribute to the process. We have 
previously noted 9] that the requirement that final states 
be outside the nuclei may not be necessary for some final 
states such as lepton pairs. The smooth cutoff conven- 
tionally implemented by a momentum dependent form 
factor in QED calculations does not seem unreasonable. 
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